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It has been found that quantum corrections can substantially affect the clas- 
sical results of tracking for trajectories close to the separatrix. Hence the 
development of a basic formalism for obtaining the quantum maps for any 
particle beam optical system is called for. To this end, it is observed that 
several aspects of quantum maps for the beam optics of spin-^ particles can 
be studied, at the level of single particle dynamics, using the proper formalism 
based on the Dirac equation. 



1 Introduction 

The theory of particle beam optics, currently used in the design and operation of 
various beam devices, from electron microscopes to accelerators, is largely based 
on classical mechanics and classical electrodynamics. Such a treatment has indeed 
been very successful in practice. Of course, whenever it is essential, quantum me- 
chanics is used in accelerator physics to understand those quantum effects which are 
prominent perturbations to the leading classical beam dynamics 0. The well-known 
examples are quantum excitations induced by synchrotron radiation in storage rings, 
the Sokolov-Ternov effect of spin polarization induced by synchrotron radiation, etc. 
Recently, attention has been drawn by Hill || to the limits placed by quantum me- 
chanics on achievable beam spot sizes in particle accelerators, and the need for the 
formulation of quantum beam optics relevant to such issues ||. In the context of 
electron microscopy scalar wave mechanics is the main tool to understand the image 
formation and its characteristics, and the spin aspects are not generally essential |4j . 

In the context of accelerator physics it should be certainly desirable to have a 
unified framework based entirely on quantum mechanics to treat the orbital, spin, 
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radiation, and every aspect of beam dynamics, since the constituents of the beams 
concerned are quantum particles. First, this should help us understand better the 
classical theory of beam dynamics. Secondly, there is already an indication that this 
is necessary too: it has been found || that quantum corrections can substantially 
affect the classical results of tracking for trajectories close to the separatrix, leading 
to the suggestion that quantum maps can be useful in finding quickly the boundaries 
of nonlinear resonances. Thus, a systematic formalism for obtaining the relevant 
quantum maps is required. This problem is addressed here for the case of spin-^ 
particle beams, at the level of single particle dynamics as the first step towards a 
more comprehensive theory. 



2 Quantization of the classical particle beam op- 
tics 

If the spin is ignored, one may consider obtaining the relevant quantum maps for any 
beam optical system by quantizing the corresponding classical treatment directly. 
The best way to do this is to use the Lie approach to classical beam dynamics, 
thoroughly developed by Dragt et at, particularly in the context of accelerator 
physics. Ignoring the effect of spin on the orbital motion, the spin motion has also 
been treated classically, independent of the orbital motion, using Lie methods . 

Let the single particle optical Hamiltonian corresponding to a classical beam 
optical system be / H(l±,P ± \ z), where z is the coordinate along the optic axis of the 
system, and r± = (x,y) and p = (p x ,p y ) represent the coordinates and conjugate 
momenta, respectively, in the transverse (x, y)-plane. We shall assume the beam 
to be moving in the positive z-direction. Then for any observable of the system, 
0{l_l,P ± ), n ot explicitly dependent on z, the z-evolution equation, or the beam 
optical equation of motion, is 

where the Lie operator : / : associated with any function of the transverse phase- 
space variables, f(rj_,p ± ), is defined through the Poisson bracket, 

. j , = y i fdf 9g df dg\ fdf dg df dg\ 
' ' 1 \dx dp x dp x dxj \dy dp y dp y dyj 

When the Hamiltonian 7i is ^-independent the solution of Eq. ([!]) can be written 
down as 

0{z f ) = exp(£: -H:)0(zi) 



2 



O ( Zi ) + £(: -H : O) ( Zi ) + (l 2 /2\) (: -H : 2 O) (*) 
+ (£ 3 /3!) (: -H : 3 O) { Zi ) + ... 

o (zi) + 1 ({-n, o}) {zi) + (£ 2 /2!) ({-«, {-n, o}}) 

+ (£ 3 /3!) ({-H, {-H, {-H, O}}}) (*) + ..., 



(3) 



relating O (zt), the value of O at an initial Zi, with 0(zf), its value at a final Zf, 
where Zf > Zi and I = (zf — Zi). When the Hamiltonian depends on z we would 
have 



0(z f ) 



P 



cxp 



dz : — Ti : 



0){zi) = (M (z f ,Zi)0) 



(4) 



where the transfer map, M. (zf, Zi), a Lie transformation, is now an ^-ordered expo- 
nential. 

To obtain the quantum mechanical formalism for the above system we can follow 
the canonical quantization rule { , } — > i[ , ] where [ , ] represents the 
commutator bracket between the corresponding quantum operators. This turns 
Eq. ([!]) into the Heisenberg equation of motion 



dO 

dz 



i r 
% 



(5) 



where the quantum Hamiltonian operator Ti,, and O for any observable, are obtained 
from their respective classical counterparts by the replacement 

( d d\ 

r ± — >r ± = r x = (x, y) , p ± — * p ± = -iHV x = I ~^^> ~ ifi d~ ) ' ^ 

followed by a symmetrization to ensure that the quantum operators are hermitian. 

From the Heisenberg picture of Eq. (||) let us go to the Schrodinger picture in 
which a wavefunction ip (r_|_; z) is associated with the transverse plane at z. The 
^-evolution of \ip(z)) is governed by the beam optical Schrodinger equation 



ih^{z)) = H\^{z)). 



(7) 



Since \ip {r±] z)\ 2 w iH represent the probability density in the transverse plane at z 
the average of any O at z will be 



(0)(z) 



dxdyip*(z)Oip(z) = {ip{z)\6\ip{z)) 



(8) 



with if) (r_|_; z) normalized as (if)(z)\if)(z)) = 1. 
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(zi)) and \ip f ) = \^(zf)), 



exp i -- 



dzH 



(9) 



The formal solution of Eq. (|7|) is, with = 

tyf) = U ( Z fi Z i) = Ufi\ipi) , U ft = p 
Thus, we get 

(6), = (6)(z f ) = (ii> f \6\il>f) = (ApfPUfM = (u}iOUfi)i ■ (io) 

From the correspondence between Eq. (|IJ) and Eq. ([|) it follows immediately that 



p 



exp 



Zf i A 

dz : — H : 



O = M (z f , z^ 6 . 



:n) 



with the definition : jJi : O 



H,0 



Note that in the classical limit, when 



: VK : O — > : — H : O, the quantum Lie transformation M. (zf,Zi) becomes the 
classical Lie transformation M. (zf,Zi). This shows that if a system corresponds 
classically to a map 



(r ±i ,P ±i ) — ► (r±f,P ±f ) = (R±(r ±i ,p ±i ),P ± (r ±i ,p ±i )) , 
then it will correspond to a map of quantum averages as given by 

(t±)i — > (L±)f = (R±(r ± ,p ± ))i , (pji — > (p ± ) f = (P ± (r ± ,p ± ))i 



(12) 



(13) 



To see what Eq. ( ]T3| ) implies let us consider, for example, a classical Lie trans- 
formation exp ( : |x 3 corresponding to a kick in the xz-plane by a thin sextupole. 
This leads to the classical phase-space map 



X j Xi , 



p f =Pi + ax- 



(14) 



as follows from Eq. This would correspond to the quantum Lie transformation 
exp(: |x 3 :) which leads, as seen from Eq. ([IB]), to the following map for the quantum 
averages: 



(x) f = (x)i , (p) f = (p)i + a(x 2 )i = (p)i + a{x)} + a((x - (x)) 2 )i 



(15) 



Now, we can consider the expectation values, such as (x) and (p), as corresponding 
to their classical values d la Ehrenfest. Then, as the above simple example shows, 
generally, the leading quantum effects on the classical beam optics can be expected 
to be due to the uncertainties in the initial conditions like the term a{(x — {x)) 2 )i in 
Eq. ([T5|) . As pointed out by Heifets and Yan, such leading quantum corrections 
involve the Planck constant % not explicitly but only through the uncertainty prin- 
ciple which controls the minimum limits for the initial conditions. This has been 
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realized earlier also P, P| [T0[|, particularly in the context of electron microscopy || |9| . 
In a detailed study of a simple example it has been found that trajectories close 
to the separatrix are strongly perturbed in spite of very small initial rms (10~ 15 ) 
and small (1500) number of turns. 

As is clear from the above, a quantum formalism derived from the classical beam 
optics can be expected to give all the leading quantum corrections to the classical 
maps. The question that arises is how to go beyond and obtain the quantum maps 
more completely starting ab initio with the quantum mechanics of the concerned 
system since such a process should lead to other quantum corrections not derivable 
simply from the quantization of the classical optical Hamiltonian. Essentially, one 
should obtain the quantum beam optical Hamiltonian 7i of Eq. (^) directly from the 
original time- dependent Schrodinger equation of the system. Once 7i is obtained 
Lie methods 0, [?| can be used to construct the quantum 2-evolution operator JJ ^ 
and study the consequent quantum maps. Derivations of 7i for the Klein-Gordon 
and Dirac particle beams will be discussed in the following sections. 

A more complete theory, even at the level of optics, must take into account 
multiparticle effects. To this end, it might be profitable to be guided by the mod- 



els developed by Fedele et al. ITT], [12] (thermal wave model - TWM) and Cufaro 
Petroni et al. [13| (stochastic collective dynamical model - SCDM) for treating the 
beam phenomenologically as a quasiclassical many-body system. Though the details 
of approach and interpretation are different, both these models suggest phenomeno- 
logical Schrodinger-like wavefunction descriptions for the collective motion of the 
beam. In TWM the beam emittance plays the role of h. In SCDM it is argued 
that h is to be replaced by an effective unit of beam emittance given in terms of the 
Compton wavelength of the beam particle and the number of particles in the beam. 
It may be noted that Lie algebraic tools can be used to handle any Schrodinger-like 
equation. 



3 Using the Klein-Gordon equation ignoring the 
spin 

One may consider getting a theory of quantum maps for spin-| particle beam optical 
system based on the Klein-Gordon equation ignoring the spin. For this, one has to 
transform the equation 



ih Wt 



+ m 1 c i ) ^ (r ± ,z;t) 



(16) 

into the beam optical form in Eq. (|7|); in Eq.([Tj]) q is the charge of the particle, 
7t_l = (ti x , 7t y ) = (p x - qA x ,p y - qA y ), ±]_ = 7r| + iri and and A = (A x , A y , A z ) 
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are, respectively, the scalar and vector potentials of the electric and magnetic fields 
of the optical system (E_ = — V0, 5 = V x 4). In the standard relativistic quan- 
tum theory |TJJ], Feshbach-Villars and Foldy-Wouthuysen techniques are used for 



reducing the Klein-Gordon equation to its nonrelativistic approximation plus the 
relativistic corrections. Applying analogous techniques in the special case of a quasi- 
paraxial (|p_J <C p z ) monoenergetic beam propagating through a system with time- 
independent fields one can reduce Eq. (|l6l) to the beam optical form of Eq. (0) with 
7i containing a leading paraxial part followed by nonparaxial parts || |9], • In this 
case the wavefunction in Eq. (Il6|) can be assumed to be of the form 



\l/(r_|_, z; t) = ip(r_j_] z) exp 



- ( Po z - Et) 



(17) 



where po is the design momentum of the beam and E = +\j c 2 p$ + m 2 c 4 . Then the 
resulting time-independent equation for ip(r ± ; z) can be regarded as describing the 
scattering of the beam particle by the system and transformed into an equation of 
the type in Eq. (Q) |, |, [15|. 

For example, for a normal magnetic quadrupole lens with A = (0, 0, ^K(x 2 —y 2 )), 
where K is nonzero inside the lens region and zero outside, 



H 



2p 



-(fx+fy 



qK [x — y" 



1 



qKh 2 



(pI-pI)- 



It must be noted that while the first three terms of 7i in Eq. (|T^) are exactly the 
terms derivable by direct quantization of the classical beam optical Hamiltonian 
the last, ^-dependent, term is a quantum correction not derivable from the classical 
theory. Though such ^-dependent terms may seem to be too small, particularly for 
high energy beams, they may become effective when there are large fluctuations in 
the initial conditions since they essentially modify the coefficients in the classical 
maps. 



4 The proper theory using the Dirac equation 

For a spin-i particle beam the proper theory should be based on the Dirac equation 
if one wants to treat all the aspects of beam optics including spin evolution and 
spin-orbit interaction. In such a case the Schrodinger equation to start with is 

ih-^9f(r L ,z-,t) = m(r L ,z-,t) , (19) 

where \1/ is now a 4-component spinor and H is the Dirac Hamiltonian 

H = f3mc 2 + q<j) + ca± ■ n_± + ca z (~^^ — qA z j — /i a /3S ■ B_ , (20) 
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including the Pauli term to take into account the anomalous magnetic moment 
fi a . In Eq. ( p0[) all the symbols have the usual meanings as in the standard Dirac 
theory |T4|]. Considering the special case of a quasiparaxial monoenergetic beam we 



can take \l/ (r ± , z; t) to be of the form in Eq. (|17D. Then the 4-component ip{z±', z ) 
satisfies the time-independent Dirac equation 



(3mc 2 + q<j) + cq_i_ ■ jt ± + ca z ( — ih— qA z ) — u a /3S • B_ 

\ dz ) 



i> (l±] z) = E^j (r ± ; z) . 

(21) 

describing the scattering of the beam particle by the system. 

Actually Eq. (^T|) has the ideal structure for our purpose since it is already lin- 
ear in So one can readily rearrange the terms in it to get the desired form of 
Eq. (|[|). However, it is difficult to work directly with such an equation since there 
are problems associated with the interpretation of the results using the traditional 
Schrodinger position operator fT(|. In the standard theory the Foldy-Wouthuysen 



(FW) transformation technique is used to reduce the Dirac Hamiltonian to a form 
suitable for direct interpretation in terms of the nonrelativistic part and a series of 
relativistic corrections. Derbenev and Kondratenko (DK) |T7| used the FW tech- 



nique to get their Hamiltonian for radiation calculations. Heinemann and Barber [|18 
have reviewed the derivation of the DK Hamiltonian and have used it to suggest 
a quantum formulation of Dirac particle beam physics, particularly for polarized 
beams, in terms of machine coordinates, observables, and the Wigner function. 

In an independent and different approach an FW-like technique has been used 
to develop a systematic formalism of Dirac particle beam optics in which the aim 
has been to expand the Dirac Hamiltonian as a series of paraxial and nonparaxial 



approximations || |9|, [TDl [19], pOfl . This leads to the reduction of the original 4- 



component Dirac spinor to an effective 2-component ip(zj_; z) which satisfies the 



accelerator optical Schrodinger equation [20 



ih^ (r ± - z) = Hi, (r ± ; z) , V *) = ( ^ % ) , (22) 

where 7i is a 2 x 2 matrix operator incorporating the Stern-Gerlach (SG) spin-orbit 
effect and the Thomas-Bargmann-Michel-Telegdi (TBMT) spin evolution. As is 
usual in accelerator theory the spin operator S_ = \TiQ_ entering the accelerator op- 
tical Hamiltonian Ti refers to the rest frame of the moving particle. Further, (x, y) 
and (p x ,Py) in IK correspond to the observed particle position and momentum com- 
ponents in the transverse plane. It should be noted that the 2-component ip{r ± ; z) of 



Eq. (|22|) is an accelerator optical approximation of the original 4-component Dirac 
spinor, valid for any value of the design momentum po from nonrelativistic to ex- 
treme relativistic region. 
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For the normal magnetic quadrupole lens the accelerator optical Hamiltonian 

reads 



(xS y + yS x ) , (23) 



where 7 = E/mc 2 and e = 2m/j la /h. The last spin-dependent term accounts for the 
SG kicks in the transverse phase-space and the TBMT spin evolution. As in the 
Klein-Gordon case of Eq. ([Tj|), Ji of Eq. (|23|) also contains all the terms derivable 
from the classical theory plus the quantum correction terms. But, it must be noted 
that the scalar quantum correction term in Eq. ( p3|) ) (4th term) is not the same 
as the 4th term in Eq. ( |l8|) . Thus, besides in the ^-dependent effects of spin on 
the orbital quantum map (e.g., the last term in Eq. (^)), even in the ^-dependent 
scalar quantum corrections the Dirac particle has its own signature different from 
that of the Klein-Gordon particle. 



5 Conclusion 

The problem of obtaining the quantum maps for phase-space transfer across particle 
beam optical systems has been reviewed. The leading quantum corrections to the 
classical maps are mainly due to the initial uncertainties and involve the Planck 
constant h not explicitly but only through the minimum limits set by the uncer- 
tainty principle. These corrections can be obtained by direct quantization of the 
Lie algebraic formalism of classical particle beam optics. The Klein-Gordon and 
Dirac theories add further subtle, ^-dependent, corrections which may become ef- 
fective when there are large fluctuations in the initial uncertainties. Contrary to the 
common expectation the scalar approximation of the Dirac theory is not completely 
equivalent to the Klein- Gordon theory. All aspects of quantum maps for spin-^ 
particle beams, including spin evolution and spin-orbit effects, can be studied, at 
the level of single particle dynamics, using the proper formalism based on the Dirac 
equation. 
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